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NEW PINE STRUCTURE COOLING RATE 


Walter R. Hoegy 
ABSTRACT 

One of the dominant electron cooling processes in the 
ionosphere is caused by electron impact induced fine struc- 
ture transitions among the ground state levels of atomic 
oxygen. Tliis fine structure cooling rate is based on theo- 
retical cross sections. Recent advances in the numerical 
cross section determinations to include polarization effects 
and more accurate representations of the atomic target re- 
sult in new lower values. These cross sections are em- 
ployed in this paper to derive a new fine structure cooling 
rate which is between 40% and 60% of the currently used 
rate. A new generalized formula is presented for the cool- 
ing rate {from which the fine structure cooling rate is de- 
rived), valid for arbitrary mass and temperature difference 
of the colliding particles and arbitrary inelastic energy 
difference. 
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INTRODUCTION 


Electron energy balance is an important aspect of a total understanding of 
the ionosphere and thermosphere. Solar EUV and energetic particle heating of 
electrons is balanced by thermal conduction and by cooling to ions and neutrals. 
Dalgarno and Degges (1968) have shown that the dominant electron cooling pro- 
cess in the lower thermosphere is due to electron impact collisions with the 
ground state of atomic oxygen which produce transitions among the spin multl- 
plets: the fine structure levels ^Pj , J = 0,1,2. A number of analytical formulas 
for this fine structure cooling rate have appeared; Dalgarno (1069), Lejeune 
and Petit (1969), and Rees and Roble (1975). Recent articles by Schunk and 
Walker (1973) and Swartz and Nisbet (1973) used the Dalgarno cooling rate 
formula. All of these cooling rates are based on the original fine structure cross 
section calculation of Brieg and Lin (1966). 

Recently, three new fine structure cross section calculations have been 
published; Saraph (1973), Tambe and Henry (1974, 1976), and LeDourneuf and 
Nesbet (1976). The importance of the fine structure cooling rate relative to 
other electron cooling rates, warrants a comparison of the cooling rates calcu- 
lated with the four independent theoretical cross sections, respectively desig- 
nated BL, S, TH, and LN for Brieg-Lin, Saraph, Tambe-Henry, and 
LeDourneuf-Nesbet, 

This paper briefly summarizes the differences in the cross section calcu- 
lations and then identifies the most accurate ones. A least square fit of the 
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theoretical cross sections is made to a simple power law for use in the cooling 
rate formula. The fine structure cooling rates are evaluated for each of the 
cross sections employing a generalized formula for the inelastic cooling rate, 
which has been derived from the inelastic Boltzmann collision operator using no 
^ approximations. Use of the general formula, Eq. 3, limits errors in the cooling 

rates to errors in the cross sections or in the approximate analytical represen- 
tation of the cross sections, Eqs. 1 and 2. The four cooling rates are compared 
in Figure 1 to the widely used cooling rate (Dalgarno, 1969), 


THEORETICAL CROSS SECTIONS 


The cross sections for fine structure transitions among the triplet J = 0,1,2 
levels of the atomic oxygen ground state have been computed independently 
by BI,, S, Til, and LN. The square of the transition matrix elements is a di- 
mensionless pure number, n, first employed by Hebb and Menzel (1940) and 
called collision strength {in analogy to line strength for radiative transitions) by 
Seaton (1953). The sirmmetry of in initial and final labels J2(J, J') = J2(J' , J) 
expresses the rule of detailed balancing, relating collisions with their inverse. 
The cross section for the transition J J' is obtained from the collision strength 
(matrix elements) by 


£7 


jj' “ 


mi“S2(J,j') 

ejkj 


cm‘ 
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where = 5.29167 x 10"^ cmisthe Bohr radius, gj = 2J + 1 the statistical weight, 
andkj theinitial wave number of the colliding electron, = 6.3335 x 10“® x E/k, 


2 



k is Boltzmann's constant and Ej the initial electron energy in ergs. The 
triplet level energies are E 2 = 0, Ej = 0. 02eV and Eq = 0. 028eV, Therefore 
the energy differences - the level splittings, when expressed in units of ®K are 
228*’K, 326”K, and98®K for the transitions 1-2, 0-2, and 0-1 respectively. A 
discussion of the four collision strength calculations follows. 

Breig and Lin made the first atomic oxygen fine structure collision strength 
calculation which led to a realization of the importance of the fine structure cool- 
ing rate. Their calculation used a single configuration for the atomic wave 
function, ls^2s^2p'^. They round that inclusion of level splitting (accounting for 
energy differences of the triplet ground state levels) increased the value of , 
Saraph used a multiconfiguration representation of the atomic wave function - 
the configurationsls^2s^2p'’ and Is22s2p‘* 3d for the ground state. In the wave 
function expansion, Saraph included the bound state ls^2s^2p^ to correctly treat 
p waves. The numerical difference between the BL and S calculations is prim- 
arily due to the absence of this important bound state in the BL calculation. All 
calculations after BL Include the bound state Is^ 2s^ 2p^. The effect of using the 
multiconfiguration representation rather than the single configuration represen- 
tation is to slightly increase the collision strengths. 

Tambe and Henry (1974) demonstrated that polarization effects (not Included 
in the S calculation) as well as level splitting decrease the collision strengths, 
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but they did not use the ls^2s2p‘’ 3d confi{?uratlon used by Saraph. The contra- 
diction with BL on the direction of the level splitting effect is explained by the 
incorrect p wave calculation of BL. Tambe and Henry (1976) used a different 
multiconfiguration representation than Saraph and demonstrated that at energies 
>3000“K, the collision strengths are not sensitive to the multi- versus single 
representation, but instead are more affected (lowered) by inclusion of a polari- 
zation potential. 

LeDourneuf and Nesbet (1976) used a more complete multi configuration 
representation (including polarization) than TH, but did not include level splitting. 
They made an approximate perturbation correction to account for level splitting 
and consequently found good agreement with the combined Tambe and Henry re- 
sults (1974 and 1976). 

The nearly identical TH and LN calculations are the best to date. They in- 
clude the most terms in the target atom representation and in the wave function 
expansions; however the low energy TH calculation should be extended to include 
a multi configuration representation, and the LN calculation should be modified 
to explicitly include level splitting. Although there are no measured fine struc- 
ture cross sections available for comparison, agreement has been found (within 
experimental scatter) between measured and calculated elastic electron-oxygen 
cross sections (Tambe and Henry, 1976). 
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COLLISION STKENGTII FITS 


Evaluation of the cooling rate requires the integral of the product of the 
cross section and the particle distribution functions over all energies. In order 
to carry this out, the collision strengths of BL, S, TH, and LN are represented 
analytically by a least square fit of the calculated points to a linear function of 
fin n versus CnT, with fit parameters A and B 

n = AT'^ (2) 

where T is the final electron energy (in ”K) in an excitation transition and is the 
initial energy in a deexcitation transition. This insures that 0 as the energy 
approaches the threshold energy (LeBourneuf and Nesbet, 1976). An alternate 
analytical fit, an expansion in T and T^ , was used by LN to represent their cal- 
culated points, however their fit gives too much weight to the high energy 
points, thus modifying some of the low energy values. Hence the fit of Eq. 2 
was used to represent all four sets of calculations. We consider Eq. 2 to be an 
approximation to a form completely symmetric in the initial and final energies 
T., Tj, where T of Eq. 2 would be replaced by \/Tj Tp. 

One problem in applying the fit of Eq. 2 to the theoretical SI values is in 
determining how the authors define their energy parameter. Representing the 
energy parameter T as the initial energy of the excitation transition instead of 
the deexcitation transition could produce differences in the cross sections of 60% 
at 500° K and 20% at 3000° K. All four calculations are consistent with the 
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assumption that initial and final energies refer to the excitation transitions. 
The fit values of A and B are listed in Table 1. 


COOLING RATE 

The electron impact fine structure cooling rate is calculated from the inte- 
gral over all energies of the product of the fine structure cross section with the 
electron and oxygen velocity distribution functions. The f irmula for this integral 
is obtained from the energy moment of the inelastic Boltarnann collision operator 
(Desloge andMatthyse, 1960). In performing this energy moment, the only assump- 
tion is that of Maxwellian velocity distributions. The result is a formula for the in- 
elastic, non-ionizing cooling rate, valid for arbitrary masses m and M and 
temperatures t and T of the colliding particles of number densities n and N 
respectively. 


di 



dxe'’“ 


■ 2k(T-t)x3/2 


•( Vx + e,j/ky (x/) 

+ (Vx + e,j/k_7 - y/x) (xy) j | 


(3) 
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where the lower integration limits are the thresholds Th, 0 for excitation, de- 
excitation transitions respectively, a is the total Inelastic cross section, is 
the inelastic momentum transfer cross section, and 

y=(niT + Mt)/(m + M) 

determines the relative concentration of the I th level, Nj , 

N, = Na, 

a, - g, Z“’ exp(-Ej/kTj) 

g, = statistical weight 

Z = partition function 

E, = level energy 

T, = level temperature, characteristic of the Boltzmann distribution of 
the levels 

Cjj = E| - Ej , energy level difference. 

Since Eq. 3 was derived with no approximations, it is applicable for any neutral- 
neutral excitation collisions, where m ~ M as well as electron impact excitation 
collisions, where For purely elastic collisions, Eq. 3 reduces to the 

standard formula for elastic cooling rates, {Banks and Kockarts, 1973), while 
for ■.^elastic collisions it contains an extra term proportional to the temperature 
difference, T-t. 

The electron cooling rate is obtained from Eq. 3 upon substitution of the 
fine structure cross section, Eqs. 1-2, neglecting higher order terms in m/M, 
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where m,t, n are replaced by the electron parameters m^, T^, n^, and M, T, N 
by the oxygen parameters M, T„, n(0), and by assuming equality of the monr.'^n- 
turn transfer and total cross sections (no differential scattering cross sections 
have been published): 


‘iUe 

cil 


= 8.629 X to-® 


n^iKQ) 

Z 


SAB! T”-'® 


• I e(D^ -E^) + 5.91 X lO'^ (T^-T,) 

• ( (1 +B)D +(E/T, + 1 + B)eJ1 , 


where the sum is over the three transitions 12, 02, 01 and in respective order 
the terms are: 


£ = 0.02,0.028,0.008 

= exp(-228/T, ), exp (-326/Tg ), cxp(-326/Tq ) 

Ej. = exp(-228/Tj.), cxp(-326/TJ, exp (" (^ + ^ )) 
E = 228, 326,98 

Z = 5 + 3 exp(-228/T, ) + cxp('-326/T(, ), 
and A and B are given in Table 1. 
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Using Eq. 4, the fine structure cooling rates for the TH and LN collision 
strengths are plotted in Figure 1 as i;atios with respect to the Dalgarno (1909) 
cooling rate 

D = 3.4X 10-'2„^n(0) (I-7X 10'® T^J. (5) 

for the condition Tg = T, = Tg = 800° K and for different values of T^. Through- 
out the temperature i-ange = 800° K to 6000° K it is seen that the LN derived 
cooling rate is ~ 10% liigher than the TH rate; tltis small difference demon- 
strates the good agreement between tlieir cross section calculations as noted by 
LeDourneuf and Nesbet (1976). The 800° K gas temperature is typical for mini- 
mum solar conditions. The ratio of cooling rates is nearly independent of Tg, 
therefore Figure 1 may be applied to a wide range of Tg values from ~ 500° K 
to ~ 1500° K. (Cooling rates based on the BL and S collision strengths are not 
shown since they did not include some important terms in their calculations.) 

Note that the fine structure cooling rate formula depends on the four tem- 
peratures T^., Tg, Tg, and Tj , It cannot always be assumed that the level temp- 
eratures Tg and Tj are equal to the gas kinetic temperature Tg. For example, 
Tohmatsu (1965) found that the observed line intensity ratios of the 1300 A day- 
glow could only be explained when Tg and Tj were significantly greater than Tg. 
To produce this temperature inequality, Tohmatsu showed that the reaction rate 
for the neutral -neutral fine structure transition would have to be orders of mag- 
nitude less than the rate estimated by Bates (1951), in such a case the electron 
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impact collisions would be sufficiently rapid to elevate and Tj above T^. 
Hence the full functional dependence of the cooling rate on all the temperatures 
has been retained in Eq. 4, Further discussion of this is deferred to a second 
paper. The term in Eq. 4 which explicitly contains T^, is of the order of magni- 
tude of the elastic electron-oxygen cooling rate. 

CONCLUSION 

The original cross section calculation by Breig and Lin (1966), although it 
neglected some important target representations, was instrumental in bringing 
fine structure cooling into electron heat balance study. Recent advances in the 
quantum mechanical collision strength calculations by S, TH, and LN have made 
possible reliable determinations of the atomic oxygen fine structure cross 
sections. The calculations of Tambe and Henry (1974 and 1976) and LeDoumeuf 
and Nesbet (1976) together represent the best calculation to date since they in- 
clude most of the effects which have been shown to be important - polarization, 
multiconfiguration representation of the target atom, and level splitting. Their 
calculations should be redone, however at energies below 3000° K to simultan- 
eously include energy level splitting and all the states considered by TH and LN, 
in order to improve the accuracy of the fine structure cooling rate. Their cal- 
culations above 3000° K are in good agreement and may be considered very re- 
liable since they are practically independent of level splitting. The dis- 
crepancy of the Saraph (1973) calculation from Tambe-Henry and LeDourneuf- 
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Nesbet can probably be attributed to the neglect of polarization effects and a 
different variational procedure used to determine the multi configuration repre- 
sentation of the target atom. 

The fine structure cooling rate has been evaluated with the new cross sec- 
tions and is given in a form, I3q. 4, applicable to the most general non- 
equilibrium condition of the ionosphere. Numerically, the new cooling rate is 
40% to 60% of the widely used Dalgarno cooling rate. Consequently, the net 
ambient electron cooling rate is reduced by almost a factor of two in the lower 
thermosphere. Thus a re-examination of previous electron heat budget studies 
is warranted by the new smaller cooling rate. Brace et al. (1976) have applied 
the new cooling rate (with level temperatures Tq and Tj equal to T^) to Atmosphere 
Explorer C satellite electron heating and cooling calculations and have demon- 
strated an over all improvement in the heat budget. 
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Table 1 

Parameters in Fit of n = AT® to Theoretical Collision Strengths 


Source 

A(12)t 

A(02) 

A(01) 

B(12) 

B(02) 

B(O.l) 

BL 

2, 531(-6)+ 

2.46{-7) 

2. 708(-7) 

1.272 

1.377 

1.414 

S 

1.561(-5) 

7.476(-6) 

3.166(-6) 

1.001 

0.9411 

1.052 

TH 

7.883(-6) 

9.466(-6) 

1. 037(-8) 

1.021 

0. 8458 

1.633 

LN 

8.58(-6) 

7.201(-6) 

2.463(-7) 

1.019 

0.8998 

1.268 


*The numbers in parentheses denote powers of ten. 
tThe numbers denote initial and final J levels 0, 1, and 2. 


14 



ONnooo sj JO oi±va 


15 


800 °K 



